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We consider and compare the structural properties of bulk TIP4P water and of a sodium chloride 
aqueous solution in TIP4P water with concentration c = 0.67 mol/kg, in the metastable supercooled 
region. In a previous paper [D. Corradini, M. Rovere and P. Gallo, J. Chem. Phys. 132, 134508 
(2010)] we found in both systems the presence of a liquid-liquid critical point (LLCP). The LLCP 
is believed to be the end point of the coexistence line between a high density liquid (HDL) and a 
low density liquid (LDL) phase of water. In the present paper we study the different features of 
water-water structure in HDL and LDL both in bulk water and in the solution. We find that the 
ions are able to modify the bulk LDL structure, rendering water-water structure more similar to the 
bulk HDL case. By the study of the hydration structure in HDL and LDL, a possible mechanism 
for the modification of the bulk LDL structure in the solution is identified in the substitution of the 
oxygen by the chloride ion in oxygen coordination shells. 

PACS numbers: 61.20.-p,64.60.My,65.20.Jk 
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I. INTRODUCTION 

The possible existence of a second critical point of wa- 
ter in the liquid supercooled metastable phase has been 
the subject of a long debate in the literature. The first 
hypothesis of its existence originated from the results of a 
computer simulation on water modeled with the ST2 po- 
tentials. On the basis of those results the thermodynamic 
anomalies of water upon supercooling were interpreted in 
terms of the long range fluctuations induced by the pres- 
ence of a second critical point. This critical point would 
be a liquid-liquid critical point (LLCP) located at the end 
of the coexistence line between a low density liquid (LDL) 
phase and a high density liquid (HDL) phase of water. 
In the LLCP scenario, these liquid phases would be the 
counterpart at higher temperature of the well-known low 
density amorphous (LDA) and high density amorphous 
(HDA) phases of glassy water. The hypothesis of a LLCP 
scenario for water motivated a large number of experi- 
mental, computational and theoretical investigations 2 ^ 3 -. 

Different interpretations of the origin of the thermo- 
dynamic anomalies of water have been also proposed as 
alternatives to the LLCP scenario. In the singularity free 
scenario^ the anomalies of water are due to local density 
fluctuations and no critical phenomena take place. Re- 
cently a critical point free scenario^ has also been pro- 
posed in which the transition between HDL and LDL 
is seen as an order-disorder transition without a critical 
point. 

A number of computer simulations, performed on su- 
percooled water with different model potentials, con- 
firmed the plausibility of the LLCP scenario^—. There 
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are also indications from experiments of the existence of 
the LLCP in bulk water— It would be approximately 
located at T ~ 220 K at P ~ 100 MPa. 

Because of the difficulties of performing experiments in 
the region where the LLCP would reside in bulk water, 
the possibility of observing the LLCP of water in aque- 
ous solutions that can be more easily supercooled^ has 
been recently explored theoretically 2 ^ and in computer 
simulation :) 22 ! 23 . Results compatible with the existence 
of a LLCP have been also found in aqueous solutions of 
salts through thermometric experiments^—. 

In a recent paper—, by means of a computer simulation 
study of the phase diagram, we indicated the possible 
detection in thermometric experiments of the LLCP in a 
NaCl(aq) solution. 

Since the detection of low and high density forms of 
water can also offer a viable path to the experimental de- 
tection of a LLCP, structural properties of supercooled 
water and aqueous solutions are of extreme interest in 
this context. The structure of water and of aqueous so- 
lutions, can be studied with neutron diffraction using iso- 
topic substitutio n 28 ' 29 or by X-ray scattering— 

In the present paper we focus on the structural prop- 
erties of bulk TIP4P water and of the NaCl(aq) solution 
with c = 0.67 mol/kg, in order to analyze and compare 
the results in HDL and LDL especially close to the LLCP. 
The paper is organized as follows. In Sec.|H]the details of 
the computer simulations are given. In Sec. IIIII we sum- 
marize the main results obtained on the thermodynamics 
of bulk water and NaCl(aq) and we present the potential 
energy of the systems. The new results for the struc- 
tural properties of the systems are presented in Sec. IIVI 
This section is divided in two parts: water- water struc- 
ture is discussed in subsection IIV A[ while the hydration 
structure of ions is addressed in subsection llVBI Finally, 
conclusions are given in Sec. [V] 
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II. COMPUTER SIMULATION DETAILS 

Molecular dynamics (MD) computer simulations were 
performed on bulk water and on NaCl(aq) with concen- 
tration c = 0.67 mol/kg. The interaction potential be- 
tween pairs of particles is given by the sum of the elec- 
trostatic and the Lennard- Jones (LJ) potentials. 
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Water molecules were modeled using the TIP4P poten- 
tial^. The details about this potential are reported in 
the Appendix. TIP4P potential is known to well describe 
water properties in its liquid supercooled stated and also 
to be able to reproduce the very complex ices phase dia- 
gram^. 

The LJ interaction parameters for the ions were taken 
from Jensen and Jorgensen^ 3 - and the ion-water inter- 
action parameters were calculated by using geometrical 
mixing rules ey = (^u^jj) 1 ^ 2 and atj = (crucjjj) 1 / 2 . The 
ion-ion and ion- water parameters are reported in Table |U 
These parameters were optimized for use with TIP4P wa- 
ter and they well reproduce structural characteristics and 
free energies of hydration of the ions^ 3 .. 

Although the presence of ions would suggest the use 
of polarizable potentials, at the moment no joint set of 
polarizable potentials for water and ions are tested as 
reliable for very low temperatures. 



TABLE I. Ion-ion and ion- water LJ interaction parameters 



Atom pair 



e(kJ/mol) 



a(A) 



Na-Na 

Na-Cl 

Cl-Cl 

Na-O 

Cl-O 



0.002 
0.079 
2.971 
0.037 
1.388 



4.070 
4.045 
4.020 
3.583 
3.561 



Periodic boundary conditions were applied. The cutoff 
radius was set at 9 A. Standard long-range corrections 
were applied for the calculation of the potential energy 
and the virial. The long-range electrostatic interactions 
were handled with the Ewald summation method. The 
integration time step was fixed at 1 fs. 

The total number of particles contained in the simula- 
tion box is Ntot = 256. For bulk water N wat = N tot — 
256 while in the case of NaCl(aq) with concentration 
c = 0.67 mol/kg, Nwat = 250 and N Na + = N cl - = 3. 
Extensive sets of simulations were run both for bulk wa- 
ter and for NaCl(aq). The range of densities investi- 
gated spans from p = 0.83g/cm 3 to p — 1.10g/cm 3 
and the range of temperatures goes from T — 350 K 
to T = 190 K. Temperature was controlled using the 
Berendsen thermostat 4 . Equilibration and production 
simulation times were progressively increased with the 



decreasing temperature. The total running times span 
from 0.15 ns for the highest temperatures to 30 ns for the 
lowest ones. The parallelized version of the DL_POLY 
package- - was employed to perform the simulations. The 
total simulation time is circa 6 single CPU years. 



III. THERMODYNAMICS IN THE 
SUPERCOOLED REGION 

In this section we start by summarizing the main re- 
sults on the thermodynamic behavior in the supercooled 
region of bulk water and of NaCl(aq) obtained in Ref. 
and then we analyze the behavior of the potential en- 
ergy of the systems studied at different temperatures. 
From extensive simulations on bulk TIP4P water and on 
NaCl(aq) with concentration c = 0.67 mol/kg, we located 
in both systems the LLCP. Details about the calculation 
of the position of the LLCP are given in Ref. [22l . 

In Fig.[T]we report a comparison of the phase diagrams 
of the supercooled regions of bulk water and NaCl(aq) as 
obtained directly from MD. Together with the LLCP of 
the two systems, we show the Widom lines, the temper- 
ature of maximum density (TMD) lines and the liquid- 
gas limit of mechanical stability (LG-LMS) lines. The 
Widom line can be considered an extension of the coexis- 
tence line in the one-phase regio n 36 ' 37 . In bulk water the 
position of the LLCP is T c = 190 K and P c = 150 MPa. 
In a very recent paper— another estimate of the LLCP 
for bulk water has been performed with a TIP4P po- 
tential with modified parameters 3 ^,. The values that the 
authors obtained for the critical point, T c = 193 K and 
P c = 135 MPa, and the Widom line appears substantially 
the same as the values that we found with the original 
TIP4P— . 

Coming back to Fig. [T] we can see that in NaCl(aq) 
the position of the LLCP moves to lower pressure and 
higher temperature, appearing at T c = 200 K and P c = 
-50 MPa. The TMD line of the solution lies circa 10 K 
below in temperature and at slightly lower pressure with 
respect to bulk water. The LG-LMS instead is almost un- 
changed with respect to the bulk. From the comparison 
of the phase diagrams for these two systems, we found 
that the main effect of the presence of the ions is to shrink 
the region of existence of the LDL22, consistent with an 
observed increased solubility of ions in HDL wate r 27 i 40 . 

Importantly, we also found that with a rigid shift in 
temperature and pressure of the phase diagram of TIP4P 
bulk water (not shown) , that brought the TMD curve to 
coincide with the experimental TMD values, the LLCP 
in bulk water appears at T c = 221 K and P c = 77 MPa, 
close to the experimental value estimated by Mishima 
and Stanley^ T c ~ 220 K and P c ~ 100 MPa. We 
note that this shift in temperature is compatible with 
the shift between the melting line of TIP4P— and the 
melting line of real water. Recently Mishima published 
a new estimate of the LLCP in the bulk at T c ~ 223 K 
and P c ~ 50 MPa 19 also compatible with our findings. 
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FIG. 1. Comparison between the thermodynamic features of 
bulk water and NaCl(aq) in the supercooled region, based 
on the results shown in Ref. I22I We report the position of 
the LLCP and of the Widom line for bulk water (dot-dashed 
line) and for NaCl(aq) (dotted line). The TMD and LG-LMS 
lines are also shown for bulk water (solid lines) and NaCl(aq) 
(dashed lines). 

These results confirm that TIP4P is a good potential for 
describing the supercooled water phase diagram. The 
same shift applied to our ionic solution led us to predict 
a LLCP in NaCl(aq) located at T c ~ 231 K and P c ~ 
— 123 MPa. These last values appear to be in a region 
accessible by experiments, being above the homogeneous 
nucleation temperature of the solutio n 20 ' 42 . 

Generally speaking, from the behavior of the potential 
energy U it is also possible to extract information on the 
thermodynamics of a system close to a phase transition. 
As already pointed out in Ref. [H and [l5| if we examine 
the curvature of the configurational part of the Helmholtz 
free energy A = U — TS, 

{dv^) T ^{dv^) T ' T {dv^l T (2) 

it must be positive in the region of stability of an homo- 
geneous phase. When the curvature of U is negative, the 
system can still be stable due to the contribution of a 
dominant entropic term in Eq. ([2]). At supercooled low 
temperatures, the stabilization induced by the entropic 
term can be less effective as the factor T in front of the 
second derivative of the entropy in Eq. [5] becomes pro- 
gressively smaller. Thus, at low temperatures the range 
of volumes where (d 2 U/dV 2 )T < 0, corresponds to a 
region of reduced stability for the homogeneous liquid, 
where separation in two distinct phases with different 
densities may occur. 

In Fig. [2] the isotherms of the potential energy are 
reported for three temperatures for bulk water and 




FIG. 2. Potential energy per molecule as a function of the 
density of the system at constant temperature. For bulk water 
(left panel) the isotherms of the potential energy are reported 
at T = 300 K, T = 245 K and T = 190 K. For NaCl(aq) (right 
panel) the isotherms of the potential energy are reported at 
T — 300 K, T = 235 K and T = 200 K. Continuous lines are 
a polynomial best fits to the simulated points. 



NaCl(aq) . We show the behavior of potential energy for 
ambient temperature T = 300 K, for the first tempera- 
ture showing a negative curvature of U, T = 245 K in the 
bulk and T — 235 K in the solution, and for the LLCP 
temperature, T — 190 K in the bulk and T = 200 K in 
the solution. In both systems, at high temperature the 
potential energy is a positively curved function of the 
density and it decreases monotonically (becomes more 
negative) as density increases. When the temperature 
is decreased, a minimum is formed after which the cur- 
vature of the potential energy becomes negative. The 
minimum corresponds to the presence of a tetrahedrally 
ordered liquid with an energetically favorable configura- 
tion^. For the temperature corresponding to that of the 
LLCP, we see that in both systems the minimum becomes 
very deep and a maximum is suggested at high densities, 
indicating the possible occurrence of a second minimum 
at higher densities. This second minimum has been pre- 
viously observed for confined bulk water— and for higher 
concentration NaCl(aq) solutions 4 ^ and connected to the 
existence of two distinct liquid phases in the system, as 
at very low temperatures the entropic contribution is de- 
pressed and the behavior of the free energy A can be 
approximated with that of the potential energy U. Com- 
paring the behavior of the potential energy of bulk water 
and NaCl(aq), we can notice that apart from the shift 
in the absolute value due the presence of ions, their be- 
havior is quite similar. Nonetheless the minimum at low 
density becomes more shallow in the solution, indicating 
that this phase is made less stable by the presence of 
ions, consistent with the fact that ions stabilize the high 
density phased 2 -. 
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FIG. 3. O-O RDFs of bulk water at T = 190 K for densities 
from p = 1.10 g/cm 3 to p = 0.86 g/cm 3 (LG-LMS). 



IV. STRUCTURAL RESULTS 

As discussed in the previous section, the study of the 
phase diagram of bulk water and of the aqueous solu- 
tion shows the presence of a LLCP in the supercooled 
region. We now discuss the structural properties of the 
systems. In the following subsection, we analyze water- 
water structure both in the bulk and in the solution. 
Then we study the hydration (ion-water) structure. 



A. Water-water structure 

In Fig. [3] we report the 0-0 radial distribution func- 
tions (RDFs) of bulk water at the LLCP temperature for 
decreasing densities from p = 1.10 g/cm 3 to the LG-LMS 
density p — 0.86 g/cm 3 . According to our phase diagram, 
the LLCP in bulk water is located at p = 1.06 g/cm 3 . As 
a general trend, there is an increase of the first peak with 
decreasing density. We note that the second peak from 
p = 0.86 g/cm 3 to p = 0.98 g/cm 3 does not substan- 
tially change position. As the LLCP is approached, from 
p = 1.02 g/cm 3 it starts to shift to lower distances. 

The LLCP is located at p = 0.99 g/cm 3 in NaCl(aq); 
thus in order to study the structural differences be- 
tween the LDL and the HDL both in bulk water and 
in NaCl(aq), we will take into account, in the follow- 
ing discussion, the RDFs at two density values which are 
well above and well below the estimated critical densities, 
namely p = 1.10 g/cm 3 for HDL and p = 0.92 g/cm 3 for 
LDL. 

In Fig. [4] we compare the water-water RDFs goo( r ), 
90H{r) and gHH(r) of bulk water obtained at the thermo- 
dynamic conditions in which water is either in the LDL or 
in the HDL region. The RDFs are plotted for T = 190 K 
at densities p — 0.92 g/cm 3 and p = 1.10 g/cm 3 represen- 
tative of LDL and HDL respectively. The RDFs at the 




r(A) 

FIG. 4. O-O (top panel), O-H (central panel) and H-H (bot- 
tom panel) RDFs for bulk water at T = 190 K. Solid lines: 
p = 1.10 g/cm 3 (HDL); dashed lines : p = 0.92 g/cm 3 (LDL); 
dot-dashed lines: p = 0.86 g/cm 3 (LG-LMS). 

density p = 0.86 g/cm 3 that represents, at this temper- 
ature, the LG-LMS of bulk water are also reported for 
comparison. We observe for the 0-0 RDF that the height 
of the first peak decreases going from LDL to HDL while 
its position does not change. The second peak instead is 
markedly different. In HDL its position shifts to lower 
distances, its height decreases and its shape broadens, 
with respect to the LDL. For the goH(r) the height of 
the first peak decreases in HDL and it moves to slightly 
lower distances. This slight shift is conserved between 
the first and the second shell, while the position and the 
height of the second shell are fairly similar for HDL and 
LDL, apart from the appearance of a shoulder between 
the second and the third shell in HDL. It can be also 
noted that while the LDL shows a very well-defined third 
shell, it disappears in HDL. For the guni?) also the first 
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FIG. 5. 0-0 (top panel), O-H (central panel) and H-H (bot- 
tom panel) RDFs for bulk water at T — 190 K and for 
NaCl(aq) at T = 200 K. Solid lines: bulk p = 1.10g/cm 3 
(HDL); dashed lines: NaCl(aq) p = 1.10g/cm 3 (HDL); long 
dashed lines: bulk p = 0.92 g/crn 3 (LDL); dot-dashed lines: 
NaCl(aq) p = 0.92 g/cm 3 (LDL). 



peak of the HDL is less intense than the one of LDL and 
slightly shifted to lower distances. It can also be noted 
a widening of the second shell in the HDL. The overall 
trend of water-water RDFs and in particular the differ- 
ence in the second shell of the 0-0 RDF clearly shows 
the disruption of hydrogen bonds between the first and 
second shell of water molecules that occurs in HDL, caus- 
ing it to have a collapsed second shell with respect to the 
tetrahedrally coordinated LDL2&. 

The results for HDL and LDL RDFs in bulk water are 
in good agreement with those found in experiments 2 ^. It 
has been shown that the TIP4P model is a very reliable 
model for the study of the thermodynamics of solid 32 
and supercooled liquid water™. Our results confirm the 



validity of the TIP4P model also for the study of the 
structural properties in the supercooled region. 

In Fig. [5] we compare the HDL and LDL water- water 
RDFs for bulk water and for NaCl(aq). In NaCl(aq) the 
overall shape and positions of the peaks remain similar 
to those of bulk water in both cases, but some differences 
can be noted for the LDL. The height of the first peak 
of goo(r), 9oH(r) and g H H{r) is slightly lower for the 
NaCl(aq) LDL with respect to the correspondent phase 
in bulk water. Furthermore in the goo{ r ) the first min- 
imum moves to lower distances in the solution and the 
height of the second peak is reduced. In the goH(f) the 
first and the second shell remain similar but the height of 
the third shell is damped in the NaCl(aq). These results 
indicate that at this concentration the effect of ions on 
water-water structure is not strong and all the features 
found in bulk water are preserved. Nonetheless, LDL 
seems to be more affected by the presence of ions than 
HDL. This can be due to a certain degree of disruption 
of hydrogen bonds induced by the ions. In fact, as shown 
in Fig. [5] we see that the effect of ions is that of reducing 
the LDL character of water- water structure and making 
it more similar to the one of the HDL. These results are 
in agreement with that observed in our study the ther- 
modynamics of these systems 2 ^. In fact, it was found 
that the range of existence of the LDL phase is reduced 
when ions are added to bulk water. 

To complete the comparison between HDL and LDL 
in bulk water and in NaCl(aq) we report in Fig. |B] the 
quantity 

h(r) = 4:Tr P r[0.092goo(r)+0A22goH(r)+0A86g HH (r)-l} 

(3) 

This correlation function has been obtained in neutron 
diffraction experiments on the amorphous HDA and LDA 
phases of wate r 45 ' 46 . The qualitative agreement shows 
the existing relation between the LDA/HDA phases of 
ice with the corresponding LDL/HDL phases of water 
as already pointed out in simulations of bulk water with 
the ST2 potential^ 5 -. We observe again that LDL is more 
influenced by the presence of ions. The modifications 
to h(r) in fact appear larger in LDL with a significant 
reduction of the peak at about 4.6 A and the different 
behavior for long distances, r > 6.5 A. This confirms the 
behavior that we observed for LDL, looking at the partial 
water- water RDFs. 

To conclude the analysis of water- water structure we 
report in Fig.[7]the oxygen-oxygen first shell coordination 
numbers as a function of density, at constant tempera- 
ture. We report them at the same temperatures that 
we show in Fig. [5] for the potential energy, T — 300 K, 
T = 245 K and T = 190 K for bulk water and T = 300 K, 
T = 235 K and T — 200 K for NaCl(aq). The first shell 
coordination numbers can be calculated integrating the 
radial distribution function until the first minimum is 
reached, and thus they represent the average number of 
first neighbors around an oxygen atom. At high temper- 
atures the coordination numbers decrease monotonically 
with the density, without reaching the value of 4 typical 
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FIG. 6. Function h(r), see text for definition. Top panel: bulk 
HDL (solid line), NaCl(aq) HDL (dashed line) and experimen- 
tal HDA (dotted line). Bottom panel: bulk LDL (solid line), 
NaCl(aq) LDL (dashed line) and experimental LDA (dotted 
line). 
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FIG. 7. O-O first shell coordination number as a function of 
density, at constant temperature, for bulk water (left panel) 
and NaCl(aq) (right panel). The temperatures reported are 
T — 300 K, T — 245 K and T = 190 K for bulk water and 
T — 300 K, T — 235 K and T = 200 K for NaCl(aq). 



of LDLi£ in the spanned range of simulated state points. 
As temperature is decreased, the isotherms of the coordi- 
nation numbers approach an asymptotic value of 4 when 
density decreases. For the lowest temperature, the co- 
ordination number reaches the value 4 already at quite 
high densities in bulk water and it actually drops slightly 
below 4 in NaCl(aq). In the case of the first shell coordi- 
nation numbers, HDL and LDL seem to be affected in a 
similar way by the presence of ions, with a slight decrease 
of the average number of first neighbors. Nonetheless, 
the weakening of the ordered tetrahedral configuration 
of LDL in the first shell induced by the ions may allow 
for a reorganization of the bonds that leads to a more 
packed, HDL-like structure, as we have seen looking at 
the RDFs. 



B. Hydration structure 

The hydration structure in NaCl(aq) at ambient tem- 
perature and in the moderately supercooled region has 
been the subject of several experimenta l 29 ' 47 ' 48 and com- 
puter simulation works— ~—. Nonetheless, the hydration 
behavior close to the LLCP has never been studied pre- 
viously. 

In order to study the general features of the hydration 
shells we show in Fig.[5]the Na- water and CI- water RDFs 
for p = 1.00 g/cm 3 at ambient temperature T = 300 K 
and in the deep supercooled region, at T = 200 K, corre- 
sponding to the LLCP temperature in the solution. For 
all the four couples we can observe that there are two 
well-defined hydration shells. The height of both the first 
and second peak tends to increase greatly upon super- 
cooling. In the case of Na-0 and Cl-0 we observe also 
for the second peak also a shift of its position to slightly 
lower distances at the lowest temperature. These results 
seem to indicate that upon supercooling the hydration 
shells tend to stabilize and to become more compact. In 
the case of the chloride ion, the Cl-0 first and second 
peaks are very well separated from the respective Cl-H 
peaks, by roughly 1 A. In the case of the sodium ion the 
distance is reduced to roughly 0.6 A and there is actually 
an overlap of the Na-0 and Na-H hydration shells. For 
the chloride ion the hydration structure of the shells is 
similar to the case of 0-0 and O-H structure (see Fig. [4] 
and Fig. [5]). Thus, we argue that chloride can substitute 
oxygen in the hydration shells forming linear hydrogen 
bonds with oxygen atoms as indicated by the sharp, first 
Cl-H peak. The Na + ion cannot simply substitute the 
oxygen atom. Its positive charge induces a major rear- 
rangement of the hydration shells that results in more 
packed structure. The behavior of the hydration shells 
of Na + and Cl~ ions is in good agreement with those 
found experimentally for ambient temperature in recent 
neutron diffraction experiment a 29 ' 47 . 

In Fig. IH1 we report the Na- water and CI- water RDFs 
for the HDL, at p — 1.10 g/cm 3 and the LDL at p = 
0.92 g/cm 3 and at T = 200 K. We can see that the hy- 
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FIG. 8. Na-water (top panel) and Cl-water (bottom panel) 
RDFs at p = 1.00 g/cm 3 and at temperatures T = 300 K and 
T = 200 K. 
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FIG. 9. Na-water (top panel) and Cl-water (bottom panel) 
RDFs for p = 1.10 fl/cm 3 (HDL) and p = 0.92 g/cm 3 (LDL), 
at T = 200 K. 




dration shells of sodium ions are not much affected by 
the HDL or LDL environment of the solvent. An in- 
crease of the first and second peak of the Na-0 RDF and 
of the first peak of the Na-H RDF is observed in LDL 
with respect to HDL but the position of the peaks re- 
mains unaltered apart from a very slight shift to higher 
distances of the second shell of Na-0 for the LDL. In 
the case of the chloride ion, the Cl-H RDF remains un- 
changed in HDL or LDL but the Cl-O RDF show some 
differences. In the LDL the first peak is higher than in 
HDL and the second peak is both higher and shifted to 
longer distances. We report in Table [Hi the first shell and 
second shell hydration numbers for the sodium and chlo- 
ride ions at the LLCP temperature. These numbers have 
been calculated by numerical integration of the first and 
second peak of the RDFs plotted in Fig. |H1 The integra- 
tion range spans form zero to the first minimum for the 
first hydration shell and from the first to the second min- 
imum for the second hydration shell. The modifications 
of the first shell in going from HDL to LDL are prac- 



tically imperceptible. The second shell appears already 
sensitive to the HDL and LDL environment as we can 
deduce from the decrease of hydration numbers in going 
from HDL to LDL. 

TABLE II. First shell, m, and second shell, ri2, hydration 
numbers for sodium and chloride ions in HDL and LDL at 
the LLCP temperature, T c = 200 K. 

Atom pair ni HDL m LDL ni HDL 712 LDL 



Na-O 
Na-H 
Cl-O 
Cl-H 



5.975 
15.912 
7.214 
6.983 



5.957 
15.364 
7.032 
7.002 



18.490 
49.474 
23.193 
33.177 



15.898 
43.521 
19.378 
27.069 



As shown in Fig. [5] and Fig. [51 the ions seem to affect 
the LDL more than the HDL. The most significant dif- 
ference in the hydration shells of the ions is noticed in 
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FIG. 10. Comparison of HDL and LDL O-O and Cl-O RDFs 
with distances rescaled by the Lennard- Jones parameter a. 



the second shell of the Cl-0 RDF. We have also already 
mentioned that the hydration structure around the chlo- 
ride ion resembles the 0-0 and 0-H structures (Fig. [5]) 
with the chloride ion able to substitute oxygen in shells 
of other oxygen atoms. To further inquire this possibility 
and to assess the effect of the chloride ion hydration on 
HDL/LDL we plot together in Fig. [10] the 0-0 and Cl-0 
RDFs for HDL and LDL with distances rescaled by the 
respective LJ interaction distance parameter a. For the 
0-0 pair a = 3.154 A as given in the TIP4P model 3 -*, 
for the Cl-0 pair a = 3.561 A as reported in Table |T| 

In both LDL and HDL the first and the second shell of 
oxygen atoms are closer in the case of the Cl-0 RDF (see 
also Fig. [S] and Fig.[S]). To make a quantitative compar- 
ison we can look at the distance between the first peak 
and the second peak of the RDFs AR. In HDL, for the 
0-0 pair AR = 0.49, for the Cl-0 AR = 0.42. In LDL, 
for the 0-0 pair AR = 0.55 and for Cl-0 AR = 0.44 
(distances in real units are obtained multiplying AR by 
the respective a). Thus we see that in LDL the chloride 
ion has the major effect in pulling inward its second hy- 
dration shell of oxygens. As shown in Fig. @] and Fig. [3] 
the main difference in the structure of the HDL and LDL 
is in the position of the second shell of the 0-0 RDFs, 
therefore, a possible explanation for why the presence of 
ions affect the LDL more than the HDL. In fact the chlo- 
ride ion can take the place of the central oxygen atom in 
oxygen shells and can be accommodated in water-water 
structure at the price of bending hydrogen bonds and 
pulling the second shell of oxygen closer to the first shell. 
This can be tolerated in HDL where the molecular struc- 
ture is already collapsed in the second shells but the 
LDL structure results instead disrupted. In essence, the 
substitution of oxygen atoms by chloride atoms together 
with the hydrogen bond disruption caused by Na + forces 
the LDL structure to become more HDL-like. Hence, we 
can understand also from a structural point of view why 



the region of existence of the LDL phase shrinks in the 
NaCl(aq) with respect to bulk water. 



V. CONCLUSIONS 

By the use of MD computer simulations, we have stud- 
ied the structural properties of TIP4P bulk water and of a 
sodium chloride solution in TIP4P water with concentra- 
tion c = 0.67 mol/kg. In particular, we were interested in 
the structural differences between the two phases of wa- 
ter, HDL and LDL, that appear in the deep supercooled 
region. From our previous work22 we knew accurately 
the phase diagrams of these systems and in particular 
the position of the LLCP. 

Our results showed that the TIP4P model can repro- 
duce well the structural properties of the two phases of 
supercooled liquid water in addition to being a very good 
potential for the study of liquid 22 and solid 32 water ther- 
modynamics. Comparing water-water RDFs in bulk wa- 
ter and in NaCl(aq) we have seen that the LDL is affected 
by the presence of ions more than the HDL, as indicated 
also by the shrinkage of the LDL phase observed in the 
study of the thermodynamics. 

The study of the hydration structure of ions in HDL 
and LDL revealed that a disturbance to the LDL struc- 
ture is induced by the substitution of oxygen by chloride 
ions in coordination shells of other oxygen atoms. The 
chloride ions in fact pull inward its second shell of oxy- 
gen atoms, disrupting the LDL structure. This, together 
with the hydrogen bond breaking caused by the sodium 
ion, causes the LDL phase to be less stable in NaCl(aq) 
solutions and its region of existence in the thermody- 
namic plane to reduce, with a consequent shift of the 
liquid-liquid coexistence line and of the LLCP to lower 
pressures with respect to bulk water. 

Since from our thermodynamic results we hypothesize 
that the LLCP region is above the nucleation line in this 
solution, an observation of the structural features pre- 
sented in this paper in X-ray and neutron scattering ex- 
periments, can also represent a viable route for experi- 
mentalists to solve the quest of the LLCP in bulk wa- 
ter. Along this line experimental indications of a HDL 
and LDL phase coexistence and LLCP have been recently 
founds. 
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APPENDIX 



The TIP4P model for water is a nonpolarizable po- 
tential where three sites are arranged according to the 
molecular geometry. The two sites representing the hy- 
drogens are positively charged with qn = 0.52 e; each 
one forms a rigid bond with the site of the oxygen at dis- 



tance 0.9752 A. The angle between the bonds is 104.52°. 
The site of the oxygen is neutral while a fourth site car- 
ries the negative charge of the oxygen qo = —2qh- This 
site is located in the same plane of the molecule at a 
distance 0.15 A from the oxygen with an angle 52.26° 
from the OH bond. The intermolecular interactions are 
represented by Eq. [1] The LJ parameters are given by 
aoo = 3.154 A and e 00 = 0.64852 kJ/mol. 
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